MIXED DARK MATTER FROM AXINO DISTRIBUTION 



S.A. Bonometto 1,2 , F. Gabbiani 3 and A. Masiero 4 



1 Dipartimento di Fisica dell'Universita di Milano, Via Celoria 16, 1-20133 Milano, Italy 
cn 2 " I.N.F.N. - Sezione di Milano 

^ 3 Dept. of Physics and Astronomy, University of Massachusetts at Amherst, Amherst, MA 01003, USA 
I.N.F.N. - Sezione di Padova, Via Marzolo 8, 1-35131 Padova, Italy 

S3 



(N 
> 
O 



Abstract 



O We study the possibility of mixed dark matter obtained through the phase space distribution of a single 
^particle. An example is offered in the context of SUSY models with a Peccei-Quinn symmetry. Axinos in 
^the 100 keV range can naturally have both thermal and non-thermal components. The latter one arises 
i^from the lightest neutralino decays and derelativizes at z ~ 10 4 . 
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I. Introduction 



Large scale structure data favour hybrid cosmological models. That dark matter (DM) consists of different 
components is an idea which dates back to the middle eighties, when the linear stages of hybrid models were 
fully evaluated®. The failure of canonical cold dark matter (CDM) models recently renewed the interest 
on such models. A recent evaluation of their non-linear stages® shows that it eases the solution of the 
problems CDM could not solve. Recent work® is based on a mixture of massive neutrinos with m u ~ 7 
eV (hot dark matter: HDM) and other non-baryonic particles yielding CDM. The HDM density parameter 
fi„ = 0.3(77^/28 h 2 eV) (h: Hubble parameter in units of 100km s -1 Mpc -1 , in this note we take h 2 = 0.5). 
Baryonic matter, expected to have a marginal role in the origin of inhomogeneities, has a density parameter 
f2e ^ 0.05, because of primeval nucleosinthesis constraints. The cold component is then assumed to 
have a density parameter I^cdm — 0.6, in order that Q = 1. Massive neutrinos derelativize at a redshift 
K) V„ 7eY). 

With respect to canonical CDM, the above CHDM hybrid model has only one extra parameter (m u ). 
Nevertheless it is hard to escape the feeling that having ^hdm/^cdm of 0(1) is a fortuitous coincidence. In 
this note we propose a hybrid model where a single particle naturally accomplishes the twofold role of CDM 
and HDM, thanks to its expected phase space distribution. The value of the ^hdm/^cdm ratio is linked to 
precise microphysical parameters. For a significant sector of the parameter space, we expect our model to 
reproduce (or to improve) CHDM outputs. It is clear that the quality of a model cannot be strictly linked 
to the request of minimizing the total parameter amount. Our choice of microphysical parameters is fairly 
natural and testable at present or forseeable experimental facilities. 

Our framework is based on the supersymmetric (SUSY) extension of the standard model where a Peccei- 
Quinn (PQ) symmetry® is implemented to solve the strong CP problem. More specifically, we consider a 
SUSY version of the invisible DFS axion®. Accordingly, beside axions, the model predicts the existence of 
its fermionic partner called axino (a). 

In order to obtain MDM, such framework is to be restricted as follows: (i) The lightest neutralino (x) 
is an (almost) pure gaugino. (ii) The SUSY soft breaking scale, related to the sfermion (/) masses is of 
O(TeV). (iii) The PQ scale (Vp Q ) is in the lower side of the interval allowed by astrophysical constraints®, 
ranging about 10 10 GeV. 

As the axino mass m 5 ~ m'j/VpQ, taking Vpq ~ 10 10 GeV and ~ TeV, leads to m a of O(100 keV). In 

turn this would be coherent only with a present CBR temperature T <C 2.7 K (or, if we require T = 2.7 K, 
with Q ^> 1, well above observational limits) unless a substantial entropy release took place between axino 
decoupling and today, e.g. at the electroweak phase transition®. Such release rises the CBR temperature 
at the observed level. 

In two previous works it was already suggested that DM consists of a's. In one of them® it was assumed 
that DM a's originated because of x decay. In the other one® it was shown that this could hardly account for 
Q, — 1, while a's, formerly in thermal equilibrium with the other components of the Universe and decoupling 
at a temperature T^dg < Vpq, can easily account for Q = 1. In this note we take into account both thermal 
and non-thermal axino components. The former component will be an effective CDM, as only fluctuations 
involving masses M ^ 0.1 M will be erased at its derelativization. The latter component, instead, can 
easily behave as HDM, derelativizing at a redshift z ~ 10 4 . Accordingly, fluctuations in such component 
will be erased up to a mass M D ~ 10 15 M . The exact values of the ratio ^hdm/^cdm and of z& er depend 
on m x (neutralino mass), m t (top quark mass), and mj. 

II. Non-thermal Axino Production 

Axinos are the natural lightest SUSY fermions; therefore can decay, according to the reaction x ~^ 5 + 7 
(we assume R-parity to be conserved). The temperature of the Universe when this occurs is given by: 

7j = Ahm 2 pl Y 2 dy /A^g x , dy . (1) 
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Here m p i = G l / 2 = 1.2xl0 19 GeV, g Xtdy — A4 + (7/8) A/"/ [A/i>(/) is the number of independent spin states 
of bosons (fermions) of particles with m < T dy , still coupled at T dy ], and 

U = ^Ln8n 3 )(myV P 2 Q )\f(y)\ 2 (2) 

is the decay rate^ 10 \ which is equated to the cosmic expansion rate. In eq. (2) f(y) = y(l — y + lny)/(l — y) 2 
and y = m 2 /m-~ (typical |/(y)| ~ 1CT 2 ). When the decay occurs, the x's have decoupled since long. This 

occurred at a temperature T X)dg when x = T Xidg /m x was already below unity: x is obtained equating the 
cosmic expansion rate with the rate of the reaction XX ~^ //■ Assuming that \ 1S a pure gaugino is crucial 
to restrict the dynamics of the reaction to pure sfermion exchange (we consider Higgs exchange contribution 
to be negligible). According to ref. (11) we obtain 

eV^-3/2 = (SalJn 2 ){Ah/2g x4g y' 2 Nmlm pl /m), (3) 

with iV = 116/27 = J2f(lf n cj (qf fermion charge; n c : color number), while g Xtdg is the analogue of g Xt dy, 
but at T dg instead of T dy . The number density of residual x's is then 



n 



0.70n rel (T x , d9 >- 3/2 e- 1/x , (4) 



where n re i(T) = (((3)/7r 2 )g n T 3 , while g n = A/5 + (3/4)A// (therefore, in our case, g n = 1.5). The number of 
a, arising from x decay, at any T <^ T Xtdy , is n d = n Xtdg [a Xtdg / a(T)} 3 , a being the scale factor, while their 
linear momentum, in average, is (m x /2)[a Xtdy /a(T)]. The detailed momentum distribution is quite far from 
thermal. It can be easily obtained taking into account that each a, when it arises (from non-relativistic 
x's), has momentum p = m x /2. Let t* be the time when a given a is produced, the comoving number of x's 
then reads 

n x = n x4g exp(-2T dy t*) . (5) 
At any time t <^ ^ dy , the momentum distribution of a's is obtainable from the relation 

drin n t* drir, dz* 

— - = -2 (6) 

dp z* dt* dp 

(redshift and scale factor are related according to 1 + z — a Q /a). Here p = m x z/2z* is the momentum of 
a's at the time t; z and z* are the redshifts at the times t and t* respectively. Owing to eqs. (5) and (6), 
the distribution of a's on momentum, at the redshift z, reads 

*,(P) = = 2-e~ q (7) 

n x4g dp p 

with q = [(z x dg / z){2p/m x )] 2 . For p = rria axinos derelativize. In average, this takes place at a redshift given 
by 

2dcr /10 4 ~ 0.2\f\(ma/keV)(m x /10 GeV)^ 2 (V PQ /10 w GeV)- 1 . (8) 

Eq. (8) immediately shows that, taking m x and m t (from which the value of / follows) in the experimentally 
admitted range, and choosing low Vpq and high mj (from which the value of m 5 follows), z dcT /10 4 ~ 1 is 
a generic consequence. Quite independently of their mass and because of their non-thermal distribution, 
axinos derelativize at the same epoch as massive neutrinos with mass ~ 7eV, and a previously thermal 
distribution. 

Together with the right derelativization, axinos also have a fair number density. Let us stress that n x 
and, therefore, n- a have a strong dependence on mj. This can be seen through eqs. (3) and (4). If powers 
of x are neglected with respect to exponentials, it turns out that n 5 oc rn~. With the large values of mj 
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considered here, x's i n the allowed mass range would close the Universe, unless they decay into a (or other 
particles with masses ^ 10 3 keV). 

Let us also notice that, once we require that m x /m a ^ 10 5 , having z der ~ 10 4 leads to z Xtdy ^ 10 9 . 
However, the value of z Xjdy is soon obtainable from eq. (2). The point is that, within the above range of 
Vpq and mj, the right value of Y dy follows from taking the top mass in the 150 GeV range, according to 
experimental constraints. 

Decay photons, emitted at z ~ 10 9 (T ~ 230 keV) are thermalised and cause no observable distortion 
of CBR spectrum. No appreciable amount of nuclides, produced in big-bang nucleosynthesis, is likely to 
be destroyed by decay photons. We shall devote the rest of this section to show this point. The essential 
argument is that, at least down to T ~ 20keV, nuclides are screened against x~decay photon disruption 
either by CBR photons or by the residual thermal electrons. Then, at T < 20 keV, the residual number of 
X is however too small to cause any appreciable consequence. 

At any given temperature T < 230 keV, the residual number density of x reads 

n x (T) ~ 10- 5 n 7 e_(230kcV/T)2 . (9) 

Here n 7 is the thermal photon number density at the same T. (If m x exceeds 30 GeV, n x would be even 
smaller.) Photons arising from x decay are well above threshold for the production of e + ,e~ pairs in the 
collision with CBR. At T ~ 100 keV the cross-section for this process is ~ 10~ 3 cr T (cry: Thomson cross 
section), while the cross-section for the reaction 7+ 2 H— * p + n is ~ 10~ 6 Or- 2 H disruption by high energy 
7's is therefore a negligible process. 

However, this is not enough to show that nuclide disruption is unlikely. In fact, in the e + , e~ pair 
production most energy goes to a single electron. Such electron, via inverse-compton in the Klein-Nishina 
regime, is going to yield almost all its energy to another photon, etc. This creates a cascade process' 12 -' and 
the low energy by-products, both electrons and 7's, are potentially dangerous to 2 H, which is no longer 
shielded by CBR 7's, as low energy by-products are mostly below threshold for e + , e~ pair production. 

As a matter of fact, at T < m e , the thermal electron abundance is depressed by a Boltzmann factor 
~ y z l 2 e~ y with y = m e /T, but, at T ~ 100 keV, a considerable amount of electrons are still present. Only 
when T ~ 20 keV (y ~ 25) does the electron number density become so low, to approach the baryon number 
density. Let us however recall that the cross section for nuclide disruption is ~ 10~ 6 <jt- The number of 
cascade by-products per %-decay is smaller than 10 6 and therefore residual electrons are a good shield to 2 H 
and other nuclides. At T ~ 20keV, the residual n x ~ 10~ 48 n r This compares with the small abundances 
of light nuclides produced in big-bang nucleosynthesis. E.g., the predicted 7 Li abundance corresponds to 
a number density ri[r Li j ~ 10~ 18 n 7 and cannot be appreciably modified by any photon cascade caused by x 
decay. 



III. Thermal Axino Production 

Aside of non-thermal a's, thermal a's, decoupling when the cosmic temperature T^^g was slightly below Vpq, 
are expected to exist' 9 ). Their present number density can be obtained from 

a Ldg T Ldg = 4-dg[^ 2 /C(^)9n,a]na,d-dg = CCfo 2 / ' C(^)9n,a\n a ,o- (10) 

(9n,a = 3/2). After a decoupling, the Universe underwent (at least) the electroweak phase transition. If we 
consider the entropy increase*- 7 -* at such transition(s), current mass limits for particles decoupling earlier are 
modified. We shall take that into account by means of a suitable p, factor in the following relation: 

Kga-d g - 9a)al_ dg Tl_ dg = g a? T^ = (^ /2)[7r 2 /C(3)]a 3 n 7 , . (11) 

Here g Q = g 1 (g l/ + <? 7 + g e )/{gj + 9e) — 43/11 (n 7jC is the present photon number density). Comparing eqs. 
(10) and (11), after simple calculations, yields the present density parameter JIcdm of thermal a. This reads 

m~ a 1 220 

ficDM = 0.51— a -— 2 - ( r. 12 

keV ph 2 [g-a-dg - 9a) 
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The spin degrees of freedom at T~ a ^ g yield g^dg — 220. Both non-thermal and thermal a's are fully non- 
relativistic today. Accordingly their present densities are p ,a-n.th — m a n x,x~dg( a x~dg/ a o) 3 and p ,a-th — 
^a n a,a-d g {.o,a-dgl a o) 3 ■ Then, owing to eqs. (4) and (12), we have that 

^HDM _ 0.70g x (gd-dg - 9a)^ x -3/2 c -l/x ^ 
^CDM QaQx-dg 

In the above relations the microphysical variables are mj, m x , mj, m t , Vpq. Other physical variables are /x, 
^cdm, ^hdm and z dcT . Besides of eqs. (8), (11), (12), the constraint rriaVpQ ~ m 2 ? also holds. 

Before concluding this section we wish to outline that the factor /x introduced here is also related to the 
baryon-to-photon ratio after the electroweak transition (where baryon number could be generated, if an 
out-of-equilibrium phase lasts long enough). In principle this might allow to relate the final value of the 
present baryon density parameter £l B to the final ratio ■ 



IV. Conclusions 

Within the frame of a SUSY implementation of the invisible DFS axion approach, we can restrict the 
parameter space either considering the expected ranges for m x , mf, m t and Vpq, or by requiring that 
z dcr ~ 10 4 , ^cdm ~ 0.6-0.8, ^cdm/^hdm ~ 3-4. The point we wish to make is that there is a significant 
overlap between the parts of the parameter space which are allowed by either set of constraints. We illustrate 
this fact by means of figs. 1-2 where microphysical parameters are constrained to provide suitable values of 
^dcr- ^cdm = 0.8 and ^hdm = 0.2 are taken, but results do not depend critically on this assumption. 

Fig. 1 shows for which values of mj and m x we obtain z de v — 10 4 . Different curves refer to the top 
masses ranging from 120 GeV to 180 GeV. Fig. 2 shows the dependence of the Vpq scale on m x , for m t = 
140 GeV and requiring z der — 10 4 . For m x ~ 30 GeV, Vpq ~ 1.7xl0 10 GeV. The corresponding value of mj 
is ~ 180 keV, while /x ~ 220. Fig. 3 is a description of the evolution of densities for different components 
in the Universe, starting from T ~ Vpq down to today's temperature. This plot is meant to show how 
complicate it can appear, a priori, to build up a component arising from heavier particle decay, with both 
the right derelativization redshift (zd er ) and final density parameter (JIhdm)- The main issue of this paper is 
that such aims are achieved fairly naturally in a SUSY implementation of the DFS invisible axion approach, 
provided large SUSY scale and small PQ scale are taken. 

The detailed evolution of density fluctuations over different scales, through a derelativization, equiv- 
alence, hydrogen recombination should be computed in order to evaluate CBR fluctuations, large scale 
structure and velocity fields, galaxy mass function. A study of the linear stages of this model will be 
presented elsewhere. It is however possible to estimate the minimal mass of fluctuations able to survive 
derelativization both for the cold and the hot parts of a distribution (Md, c and Mo,h, respectively), just by 
evaluating the mass scales entering the horizon at the time t of derelativization. 

For the cold component, M D:C ~ (27r 3 /45)(?aT|t 3 , while t = (45/87r 3 5f ) 1 / 2 m pi T -2 ; the relation between 
a temperature and photon temperature is T a = (g a / '^_ (ig /x) 1/ ' 3 T '. Accordingly, requiring T a = m a , we have 
that 

M DC ~f^y /2 ^^& d4) 

lvl D,c — \ 1O0 _ 3 j „ 2 (j2 m 2 I 1 V 



128tt 3 7 gi_ dgl x 



is ~ O.IMq. Such small value of Md )C is due to the actual a temperature. 

On the contrary M D>h ~ (tt 2 / 4Q)t z^J 2 ~ 3.6 x lO 15 (lO 4 /2;dcr) 1 ' 5 M0. Non thermal a distribution has 
a wide energy spread. Derelativization can be therefore expected to occur more gradually than for standard 
7eV neutrinos. 

In the context of a phenomenologically viable SUSY model it is therefore possible to obtain a phase 
space distribution of a single particle which provides a valid hybrid cosmological model. 
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Figure Captions 

Fig. 1 - Values of mf and m x leading to z der ~ 10 4 for f2 CDM = 4f2 H DM = 0.8. Different curves (bottom to 
top) refer to m t from 120 GeV to 180 GeV. 

Fig. 2 - Vpq as a function of m x for m t = 140 GeV and z der ~ 10 4 . 

Fig. 3 - Evolution of densities of radiation (solid line), CDM (dashed line) and HDM (dotted line) in 
logarithmic scales. The HDM line accounts for x's prior to their decay and then for non-thermal 
axinos. The large increase of pT~ 4 at T ~ 100 GeV is related to entropy release at the electroweak 
phase transition. 
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